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ABSTRACT 

In protostellar cores where the dust temperature has been raised above 100K and 
subsequently allowed to fall below the condensation temperature of methanol, recon- 
densation on to cooling grains removes methanol molecules from the gas at rates faster 
(~ 1000 times) than those of chemical reactions. Molecular recondensation can have a 
profound effect on the chemical composition of hot cores. The methanol chemistry of 
hot cores is solved analytically and the trend in molecule binding energies required is 
compared with theoretical and experimental values. 

It is demonstrated, through a model calculation incorporating recondensation, that 
it yields a consistent explanation of the similarity of the CH 3 OH , CH 3 OCH 3 and 
HCOOCH 3 abundances measured in G34.3 and W3(H 2 0) by Macdonald et al. (1996) 
and by Helmich & van Dishoeck (1997). These observations suggest that the latter 
molecules could not be derived from CH 3 OH through gas phase reactions. The manner 
in which molecular recondensation could affect the interpretation of hot core chemistry 
in general, particularly of organic molecules, is briefly discussed. 


Subject headings: ISM: clouds - ISM: abundances - ISM: molecules - molecular processes 
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1. Introduction 

In the energetic environments of forming massive 
and low-mass protostars, thermal heating or sput- 
tering in shock waves can liberate the molecular 
ice mantles that cover interstellar dust grains (e.g. 
Blake et al. 1987; van Dishoeck et al. 1993). Ob- 
servations of these hot cores can therefore be used 
to to probe the nature of grain surface chemistry. 
The desorbed molecules can themselves initiate dis- 
tinctive gas phase reaction pathways leading to the 
growth of highly complex organic molecules (Charn- 
ley et al. 1995), oxygen-nitrogen molecular differen- 
tiation (Charnley et al. 1992; Caselli et al. 1993), 
and the possibility of comparing dynamical time- 
scales through Molecular clocks’ (Charnley 1997a; 
Hatchell et al. 1998). One success of this scenario 
is the explanation, based on reactions originally pro- 
posed by Blake et al., of the high observed abun- 
dances of dimethyl ether (CH 3 OCH 3 ) and methyl 
formate (HCOOCH 3 ) as being due to protonation 
of methanol molecules evaporated from grain sur- 
faces followed by methyl cation transfer reactions with 
methanol and formaldehyde leading, following elec- 
tron dissociative recombinations, to CH 3 OCH 3 and 
HCOOCH 3 respectively. 

Recent line surveys of star-forming regions have, 
however, raised problems for the standard theory 
of hot core chemistry (Brown et al. 1988; Charn- 
ley et al. 1992). Observations of the G34.3 and 
W3(HoO) cores indicate that the CH 3 OH abundances 
are lower than in other cores and comparable with 
those of its putative daughter molecules (Macdonald 
et al. 1996; Helmich & van Dishoeck 1997). Table 
1 summarises the relevant observational data. As 
models predict that at most a few per cent of the 
evaporated methanol is converted to CH 3 OCH 3 and 
HCOOCH 3 (Charnley et al. 1992), it would be diffi- 
cult to reproduce the CH 3 OCH 3 abundance in G34.3 
even with the higher CH 3 OH abundance determined 
by Mehringer & Snyder (1996). It has therefore been 
suggested that they cannot be formed by gas phase re- 
actions, as previously thought, but that they are, like 
methanol, the products of grain surface reactions. 

While the proposed gas phase synthesis of HCOOCH 3 
is not supported by methyl cation transfer experi- 
ments (see Karpas & Mautner 1989), and in fact may 
involve HCOOH instead of H 2 CO (Charnley 1999a), 
the situation is much more serious for the currently 
accepted origin of CH 3 OCH 3 . The reaction 


CH 3 OH+ + CH 3 OH — . CH 3 OCH 4 + H 2 0 (1) 

is well-studied in the laboratory with a measured 
rate coefficient (Karpas ic Mautner 1989) and its 
inclusion in models of methanol-rich hot cores re- 
produces the CH 3 OCH 3 /CH 3 OH abundance ratios 
typically observed in many sources (Charnley et al. 
1992; Caselli et al. 1993). It is therefore puzzling 
that this theory should be at variance with the ob- 
servational data for the G34.3 and W3(H 2 0) cores. 
There is apparently insufficient CH 3 OH to produce 
the CH 3 OCH 3 and HCOOCH 3 by ion-molecule chem- 
istry, suggesting that the standard theory is in need 
of revision and perhaps indicating a mantle origin for 
these molecules. 

A key assumption of all previous hot core models 
has been that once the dust grain temperature, 7d , is 
raised above the condensation temperature of water, 
and the entire mantles evaporated, the grains remain 
hot, with 7d > 100K, for the nominal ‘age’ of the core, 
estimated to be ~ 10 5 years (e.g. Walmsley 1989). 
Theoretically, the gas and dust in hot cores should 
be thermally well-coupled with T ~2d ~ 100 -300K; 
observations by Schilke et al. (1992) show that this 
is indeed the case in Orion-KL. In these dense re- 
gions (mr — 10 6 — 10 8 cm“ 3 ), the gas-grain collision 
rates for neutral molecules X , A(X), are high. How- 
ever the associated thermal desorption rates, £(X), 
are very large for 2 d above 100 K since they have an 
exponential dependence on T<j 

£(X) = i/ x exp[-^-] (2) 

where ux is the vibrational frequency of X in a 
surface binding site, typically 10 12 - lO^s” 1 , k is 
the Boltzmann constant, and E\ is the binding (or 
adsorption) energy for phvsisorption (e.g. Sandford 
& Allamandola 1993). This is the reason why hot 
core chemistry is highly transient and requires a time- 
dependent treatment. For cooler cores with 2d below 
about 40K, the £(X) for the most abundant volatiles 
(i.e. X= CO. 0 2 , No ) are sufficiently low that a 
gas-grain chemical steady-state can occur (Charnley 
1997b). Hence, in a core where the dust has been 
raised above 100K and subsequently allowed to cool, 
molecule X will condense back on to the cooling grains 
when 2 d falls below its condensation temperature. 

In this paper I consider the effect of molecular 
recondensation on hot core chemistry. I describe 
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the relevant gas phase chemistry of CH3OH and 
CH3OCH3 semi-analytically and show that although 
methanol molecules can be removed from the gas at 
rates comparable to, or in excess of (~ 1000 times), 
those of chemical reactions, CH3OCH3 molecules can 
still be present in the gas. The detection of solid 
HCOOH by the Infrared Space Observatory (Schutte 
et al. 1996) at about 1% of H2O would correspond to 
a gas phase abundance upon evaporation of ~ 10“ 7 . 
However, this is in conflict with gas phase HCOOH 
abundances actually observed in star-forming regions 

10“ 10 — 10“ 9 ) and so I also explore whether re- 
condensation of HCOOH could resolve this issue. 

2. Methanol and Dimethyl Ether in Hot Cores 

In this section I first present the relevant details 
of the methanol (MeOH) and dimethyl ether (Me 2 0) 
gas phase chemistries, showing quantitatively how the 
Me20/MeOH ratio could be affected by recondensa- 
tion. Once evaporated, methanol is protonated by 
H 3 0 + , HCO+ and H3 but is reformed by dissociative 
recombination with electrons. It is removed primarily 
through reaction (1) and by He + and C + . Utilising 
the reaction rate coefficients compiled by Millar et al. 
(1997a), the MeOH-Me20 gas phase chemistry obeys 
the following differential equations for the fractional 
abundances y(X) = n(X)/nn (with n(X) the number 
density of X) 

nJJ l y(MeOH) = 300T 3 “ o 5 x e y(MeOH+) 

- 0.076T3- 1 6 y(MeOH)y(MeOHt) 

+ [2.5y(H 3 0 + ) 4 2.7y(HCO+) 4 5.7y(H+) 

+ 2.2y(He + ) + 7.3y(C + )]y(Me0H) (3) 


n^ l y( Me 2 0) = ISOTa’ 1 6 x,y(Me 2 OH+) 

- {[2.7y(H 3 0+) + 2.1y(HCO + ) 4- 2.0y(H+)] 

4 2[y(H+) 4 y(He + ) 4 y(C + )]}y(Me 2 0) (4) 

where y(X) = dy(X)/cft & , and dt$ = 10“ 9 df. The 
electron fraction has been written distinctly as x e and 
T 3 = T/300K. After about one year post-evaporation 
the gas ionization is in a steady-state to a good ap- 
proximation. When no metals (e.g. Na, Mg, Fe) are 
evaporated the electron fraction is given by 


x e % y(Me0H+)4y(HC0+)4y(H 3 0 + )4y(H+) 


and is constant for the age of the core. The distri- 
bution of positive charge is determined by a proton 
transfer cascade from H 2 to CO to H 2 0 , terminat- 
ing at MeOH when NH 3 /MeOH< 1. As H 2 0 and 
MeOH are removed, HCO + becomes the major ion at 
longer times (cf. Chamley et al. 1992,1995). The loss 
of protonated methanol is governed by ion-molecule 
reactions rather than grain collisions and its steady- 
state abundance can be expressed as 


y(MeOH+) _ 2.5y(H 3 0 + ) 4- 2.7y(HCO+) + 5.7y(H+) 
y(MeOH) 300T 3 -° 5 x e + 0.076r 3 -‘ 6 y(MeOH) 

( 5 ) 

Substituting equation (5) in equation (3), we can 
write the differential equation for the methanol abun- 
dance due solely to gas phase chemistry, it is 


y(MeOH) = ([^- 


a — y(MeOH) 
y(MeOH) 


- 1 


6 — cjy(MeOH) 


where 


( 6 ) 


a = 3.95 x 10 3 r 3 M * e 


b = n H [2.5y(H 3 0 + )4 2.7y(HC0 + )4 3.1y(H+)] 
c = n H [2.6y(H+)4 2.2y(He + )4 7.3y(C + )] 


The appearance of y(Hj ) in both b and c is be- 
cause H3 protonates MeOH both dissociatively and 
nondissociatively. For standard values of the cos- 
mic ray ionization rate, C ~ 10“ 17 s” 1 and a core 
density of hh ~ 10 7 cm” 3 , the electron fraction is 
x e ~ 10“ 9 (Charnley et al. 1995) and so at 200K. 
c ^ 2.5 x 10“ 6 . The maximum values of y(MeOH) 
obtained from both gas phase and solid phase obser- 
vations are around 5 x 10“ 7 (Menten et al. 1988: 
Schutte 1996) and so equation (6) has the solution 
y(MeOH) oc exp [— (564 c)^] early in the evolution, 
and the solution y(MeOH) oc exp[— cfg] at longer 
times when a > y(MeOH). 

We now consider the effect of the gas-grain inter- 
action. The collision rate of methanol molecules with 
the dust grains is 

A(MeOH) = 4.55 x 10- 18 f— V ’;i H (7) 

' MeOH ' 


where MmcOH is the molecular weight. In writing 
equation (7), unit sticking efficiency on 0.1pm grains 
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and a dust number density of 10 “ 12 hh have been as- 
sumed. Using the same scaling as for the chemical 
rate coefficients in equations (3) and (4), one has 

A 9 (MeOH) = 1.4 x KT 8 ^ 0 5 n H (8) 

The relative importance for methanol loss through 
chemical reactions as opposed to sticking collisions 
with dust grains is given by the ratio 


chemistry has evolved to some post-evaporation time 
t such that an appreciable abundance of Me 2 0 has 
been formed via reaction (1) and the methanol abun- 
dance is still close to that originally injected. If all 
the methanol molecules were to be instantaneously re- 
moved, the Me 20 abundance in the gas would not be 
affected instantaneously, but instead Me 20 molecules 
would be removed after r seconds where 


c _ 2.6y(H+) + 2.2 y(He+) + 7.3y(C+) 
A 9 (MeOH) 1.4 x 10- 8 T 3 ° 5 

(9) 

Noting that, for these ions, loss in gas phase reactions 
dominates that by grain collisions, we can evaluate 
equation (9) using analytical formulae derived from 
Charnley (1997b). 


n(He + ) 


0.5C 9 y(He) flQ v 

1.6y(CO) + 1.8 x 10” 5 exp(-35/T) 


n(C + ) « 


1.6y(He + )n(C0) 

2.0y(H 2 O) 




0.5<s 

1.7y(CO) 


( 11 ) 

( 12 ) 


where is 10 9 C- Adopting a core density of mh = 
2 x 10 7 cm“ 3 , T=200K,C9 = KTV 1 , y( CO) = 1(T 4 , 
y(He) = 0.1 and y(H 2 0) = 10~ 5 , one finds that 
c/A 9 (MeOH) » 6x 1CT 4 . Hence, if the dust grains are 
at a temperature below the condensation temperature 
of methanol, almost all colliding methanol molecules 
are retained on the grain surfaces - recondensation 
of MeOH dominates the removal of methanol over 
that due to chemical reactions. This adds a term 
-nJJ l A 9 (MeOH)y(MeOH) to the right-hand side of 
each of equations (3) and (4). The solution of equa- 
tion (6) becomes y(MeOH) a exp[-(c + A 9 )f 9 ] « 
exp[— Xt] 

What effect will MeOH recondensation have on the 
evolution of its daughter molecules? If the binding en- 
ergy of MeOH on dust is much greater than that of 
Me 2 0 and the grains are sufficiently warm that Me 2 0 
does not recondense, we can answer this question by 
considering the Me 2 0 differential equation (4). The 
production of Me 2 0 is intimately connected to the 
abundance of MeOH through equation (5); it is also 
destroyed by the low-abundance ions H + , He + and 
C + but, unlike MeOH, Me 2 0 is only produced in half 
of the electron dissociative recombinations of t he pro- 
tonated form, Me 2 OH + . Now, suppose that the core 


r = — {2.7y(H 3 0+) + 2.1y(HC0 + ) + 2.0y(H+) 

+ 2[y(H + ) + y(He + ) + y(C + )]} -1 (13) 

When methanol is absent, H 3 0 + is the dominant 
ion, y(H 3 0 + ) ft x e % 10“ 9 , and so, for a core den- 
sity of hh ~ 10 7 cm“ 3 , one finds r(Me 2 0) « 10 3 
years. Hence, there should exist values of T& and t 
for which, as the grains cool, MeOH is selectively re- 
condensed but the MeOH already incorporated into 
Me 2 0 can remain in the gas at a reasonably high 
abundance, until eventually it too is removed, ei- 
ther by recondensation, if the grains cool sufficiently 
rapidly, or by chemical reactions. 

3. Dust Temperature Fluctuations Around 
Massive protostars 

Molecular recondensation, as outlined above, clearly 
requires time variation in the local dust grain temper- 
ature within a hot core. Grains must be heated suffi- 
ciently to desorb H 2 0 , CH 3 OH , and also HCOOH, 
then be allowed to cool to about 90-100K. Luminos- 
ity variations, L m (t), of the central protostar offer the 
simplest means of realising this. A dust grain located 
a distance r from the protostar, where the emission is 
optically thin, will have a temperature (e.g. Emerson 
1988) 


T<i{r,t) oc (14) 

where the flux-averaged emissivity is assumed to have 
a v 1 frequency dependence. The luminosity of a mas- 
sive protostar can be written as (cf. Stahler, Palla k 
Ho 1999) 

L.(t) = + Lin c (15) 

H* 

where M m and R m are the protostellar mass and 
radius, .U ( t ) is the mass accretion rate and L in t is 



the intrinsic luminosity due to protostellar contrac- 
tion at the end of the convective phase. As Lint 
ocAf* R * , this term comes to dominate L m (t) as 

. the mass of the protostar increases! However, during 
the earliest phases of formation, when A/* is small, 
L m (t) is set mostly by the accretion luminosity. Tem- 
poral variations in A/(t) could therefore potentially 
have the most pronounced effect on L.(t), and hence 
on 7d(r, 0« Mass accretion must start to diminish at 
some point. Stahleret al. (1999) suggested that M(t) 
may decay rapidly relative to changes (i.e. growth) 
in Li n t. The mechanism responsible for the ultimate 
cessation in M(t) is uncertain but could be linked to 
the presence of outflows or to magnetic support of the 
parent cloud core (see Stahler et al. 1999). 

One may therefore speculate that, as a massive 
protostar begins to form, sporadic variations in M{i) 
could occur throughout the earliest phases of hot core 
evolution, mediated by the same processes that ulti- 
mately terminate mass accretion. A more detailed 
discussion of this point is beyond the scope of this 
paper, however, we note that, although having differ- 
ent physical origin, strong variations in both M(t) 
and L m (t) are known be important in the earliest 
phases of low-mass star formation - the FU Orionis 
phenomenon (Hartmann et al. 1993). As massive 
protostars spend much of their youth still embedded 
in their natal core, small fluctuations in their lumi- 
nosity are much more difficult to discern. As only a 
lowering of T d by a few kelvin may be sufficient to 
permit recondensation, only a modest reduction in 
M(t) is necessary. For example, a reduction in M(t) 
by about a factor of 5 will lower r d (r,f) by a factor 
of 1.5. 

4. Sui’face Binding (Adsorption) Energies 

We now address the question of what absolute and 
relative E\ values for each of these four molecules 
will make the recondensation scenario tenable, how 
consistent these are with known experimental or the- 
oretical values, and the associated implications for the 
surface structure of interstellar grains. We adopt the 
convention of expressing binding energies in temper- 
ature units. 

The preceding analysis considered only gas phase 
chemistry and, as methanol removal cannot occur in- 
stantaneously in practice, one must consider molec- 
ular desorption from grains. The methanol gas- 
grain chemistry is governed by the coupled differential 


equations 

y(MeOH) = n H G(MeOH) - A 9 (MeOH)y(MeOH) 
+ ^ex P [-^gp]y(MeOH) (i 6 ) 


y(MeOH) = A 9 (MeOH)y(MeOH) 

- ^exp[-^p]y(MeOH) (17) 

where g( MeOH) is the fraction of MeOH molecules 
resident on grain surfaces, G(MeOH) is the right-hand 
side of equation (3) and v 9 = 10 9 z/ M eOH =2.2xl0 21 s _1 
(Sandford ic Allamandola 1993). The abundances of 
HoO and Me20 and HCOOH obey similar equations. 
The exact chemical composition at any T& will depend 
crucially on the relative values of the binding energies, 
E\ adopted. 

The recondensation theory implies several relations 
must hold. It is intuitively obvious that one requires 
£m«oh > Em* 2 o- In both W3(HoO) and G34.3 water 
evaporation has occurred, and each still exhibits high 
gas phase abundances (Table 1; Gensheimer et al. 
1996) indicating that substantial H 2 O recondensation 
has not occurred in either core. These observations 
therefore require #MeOH > Eh 2 o Compared to their 
maximum (i.e. grain mantle) abundances of ~ 10" 7 
and ~ 10“ 6 , formic acid appears to be more depleted 
than methanol in hot cores where the methanol abun- 
dance is only slightly less than its mantle value, prob- 
ably due to destruction by gas phase reactions. If the 
inferred HCOOH depletion is in fact due to reconden- 
sation then we require L’hcooh > £ M eOH- In tota, » 
the theory suggests the following hierarchy of binding 
energies 

^HCOOH > L^MeOH > £h 2 o » £m C3 o (18) 

It should be emphasised that the conditions under 
which recondensation may occur in hot cores are dif- 
ferent from those whereby gas species accrete on to 
dust in dark clouds. Apart from the higher gas and 
dust temperatures, the probable nature of the surface 
for sticking needs to be considered. 

Early in the evolution of the cold gas that ulti- 
mately is incorporated in a hot core, molecular species 
accrete onto refractory grain cores of silicate, graphite 
and amorphous carbon, or congloinorations of these. 
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At 10 K, atoms can be physisorbed and so a mantle of 
primarily H 2 0 molecules is first to form (e.g. Brown 
k Charnley 1990). The appropriate binding energies 
are those of molecules sticking to 10K surfaces of ei- 
ther mainly H 2 0 , or, to CO, 0 2 and N 2 if a volatile 
crust surrounds the water molecules. In fact, astro- 
chemical gas-grain models of cool, ultradense cores 
(Charnley 1997b) and protostellar disks (Charnley 
1995; Aikawa et al. 1997) show that the binding 
energies needed to maintain a viable gas phase are 
close to those determined for physisorption between 
pure nonpolar or weakly-polar substances, e.g. CO 
on CO, consistent with the outermost layers compris- 
ing mostly of these molecules and not much H 2 0 . As 
the core heats up, reaching in excess of 100K, out- 
gassing is controlled by these binding energies, i.e. 
those measured by Sandford k Allamandola (1993), 
and molecules are gradually evaporated according to 
their relative binding energies. 

In hot cores, recondensation may occur at dust 
temperatures of about 90-100K and the appropriate 
surface for the initial sticking collisions will, accord- 
ing to (18), not be a water ice mantle. There are 
two possibilities: either bare silicate and/or graphite- 
amorphous carbon cores, or the organic residue sur- 
rounding these refractory cores after millions of years 
of ultraviolet photolysis and radiolysis in the molec- 
ular cloud (e.g. Allamandola et al. 1997; Moore k 
Hudson 1998). Hence, the values of Ex applicable 
may be quite dissimilar from those usually considered 
for cold molecular clouds. 

The published data falls roughly into three cate- 
gories. First, condensation temperatures derived from 
vapour pressure data for enthalpies for pure molec- 
ular gas undergoing transitions from the liquid or 
solid phase (Nagagawa 1980; Yamamoto, Nagagawa k 
Fukui 1983; Leger, Jura k Omont 1985). Physically, 
the adsorption enthalpy in this case should corre- 
spond to the surface binding energy for that on a pure 
crystal of like molecules. Second there are theoreti- 
cal estimates of surface binding/adsoprtion energies 
for molecules X on a silicate surface (Allen k Robin- 
son 1977; Hasegawa k Herbst 1993). Finally, there 
are direct measurements of Ex from temperature- 
programmed ice desorption experiments (e.g. Sand- 
ford k Allamandola 1993). 

Physisorption involves weak van der VVaals-London 
bonding and E\ depends on the polarizability of the 
adsorbed species, q\, and the composition of the 
surface, S, through the polarizability of the surface 


atoms. One finds Fx aa xcts (Masel 1996). Table 2 
summarises the relevant theoretical and experimen- 
tal data on binding energies, as well as some relevant 
molecular polarizabilities. It is clear that the theoret- 
ical binding energies do not reproduce the desired in- 
equalities of (18) and in fact would predict that Me 2 0 
should be the most susceptible to recondensation. 

These theoretical values are also at variance with 
the experimental results. The Ex for adsorption 
on silicates, listed by Allen k Robinson (1977) and 
Hasegawa k Herbst (1993), are extrapolated from 
studies of hydrogen adsorption graphite. London’s 
theory is succesful in predicting physisorption ener- 
gies for nonpolar molecules on nonmetallic surfaces, 
such as C0 2 on graphite, but extensions of the theory 
to polar molecules are generally less successful (Masel 
1996). Table 2 shows that for C0 2 the theoretical and 
experimental values are in accord. However, this is 
not the case for H 2 0 adsorption where there is a large 
discrepancy. Such a low Eh 7 o could not reproduce 
the H 2 0 ice column density observed in the envelopes 
of OH/IR stars, environments where water molecules 
are actually condensing on to warm, bare, silicate 
grains; in this case one requires Eh 2 o > 4000K 
(Meyer et al. 1998). 

The presence of the -OH group on each of the 
molecules of interest point to the importance of hy- 
drogen bonding in determining the JFx- The strengths 
of hydrogen bonds are typically 10-40 kJ mol" 1 , com- 
pared to a few kJ mol -1 for van der Waals bonding. 
This is the reason why the experimental values for 
water and methanol are so much larger than the the- 
oretical ones (cf. Table 2). The importance of hy- 
drogen bonding in this context was noted by Leger et 
al. (1985). Neglect of hydrogen bonding is the reason 
why the theoretical binding energies closely track the 
molecular polarizabilities. 

For hydrogen bonding of HCOOH and MeOH to 
occur, other OH-containing molecules must already 
be present on the grain surface. There are two cases 
of interest. First, bare silicate grains may have a sur- 
face structure similar to silanol with regularly spaced 
Si-OH groups available for hydrogen bonding (Tie- 
lens k Allamandola 1987). As the surface will be 
damaged, the occurrence of OH groups will be ran- 
dom in practice, and some of the Si atoms could also 
be replaced by others in the solid (e.g. Fe and Mg). 
The actual adsorption potential at any point on the 
surface may depend on t he proximity of a hydroxyl 
group with a contribution from the polarizability of 



adjacent species. The relative polarizabilities are in- 
deed in the desired ordering (cf. equation 18) 

<*HCOOH > CkMeOH > <*H 2 0 (19) 

Second, an organic residue containing OH groups 
could be present. The ices in hot cores are rich 
in aliphatic alcohols prior to evaporation (Charnley 
et al. 1995; Charnley 1999a). If high molecular 
weight alcohols are not desorbed significantly when 
7d > 100K, then they might provide the necessary 
binding sites. Apart from hydrogen bonding amongst 
themselves, these compounds have larger polarizabil- 
ities and so are not easily sublimed. Support for this 
picture comes from the experiments of Bernstein et 
al. (1995) who found that EtOH was retained at 150- 
200K, after H 2 O and CH 3 OH had evaporated, in the 
organic residue that persisted at the end-point of UV- 
photolysis and heating. This implies £ , EtOH>£ ? McOH 
and suggests that the low ethanol abundance seen in 
most hot cores could also be due to recondensation. 
The random occurrence of surface hydroxyls and the 
orx distribution probably also influence the surface 
potential, as remarked above. It is perhaps worth not- 
ing that the strength of the hydrogen bond amongst 
carboxylic acids is stronger than amongst alcohols, 
and both are stronger than water (Pauling 1940). If 
this trend carries over for binding to arbitrary surface 
OH groups, then it would also reproduce the desired 
E x ordering. 

5. A Numerical Model 

This effect can easily be demonstrated in hot core 
models and here a calculation is presented to illus- 
trate the ideas of the preceding sections. For con- 
stant ion fractions, equation ( 6 ) is an Abel equation 
of the second kind and although it is possible to solve 
many aspects of the chemistry analytically (Charn- 
ley 1999b), this is generally not possible when T d is 
a function of time, T^t). The results presented here 
were therefore obtained using a numerical model. 

5.1. Astronomical Binding Energies 

We need to estimate what numerical values of Ex 
are inferred by the observations. We assume that all 
the H 2 0 , CH 3 OH and HCOOH present in either gas 
of ice originated in grain surface chemistry. The total 
abundance of molecule X in a core is then 

Y(X) = y(X) + g(X) (20) 


Unfortunately the ice inventory towards each of the 
three sources in Table 1 is not known but we can take 
‘canonical* values for the Y(X) as measured in other 
sources. Relative to molecular hydrogen (i.e. 2y(X) 
etc.), for H 2 0 , MeOH, and HCOOH these are 5 x 
10 ~ 5 , 1 x 10~ 6 and 2 x 10 “ 7 , respectively (Gensheimer 
et al. 1996; Blake et al. 1987; Menten et al. 1988; 
Schutte et al. 1996). Ignoring gas phase chemistry, 
we can balance accretion and desorption rates to get 

* = r.hpjxw , '(f5) u S(^_i) 1 

( 21 ) 

The actual dust temperatures are also not accu- 
rately known but T sets an upper limit. Adopt- 
ing 7d=100K, using from Sandford ic Allaman- 
dola (1993), and the values for G34.3 in Table 1 , we 
find £ , h 3 o=517T K. Using this binding energy and 
the parameters of W3(H 2 0) gives T d = 95 K, again in 
good agreement, and so we adopt the above value for 
£h 3 o- This is in excellent agreement with the mea- 
sured value for annealed H 2 0 (Table 2) and is consis- 
tent with the grains being gradually heated and an- 
nealed prior to evaporation with any recondensation 
yielding crystalline ice. Indeed, the 3/rm water ice 
absorption spectrum of the BN object in Orion does 
show evidence for a crystalline component (Smith et 
al. 1989). 

Assuming 95K dust in W3(H 2 0) leads to £ M eOH=5264 
K. In G34.3, with 100K dust, Em<sOH is estimated 
to be in the range 5295-5410 K, depending on which 
methanol abundance in Table 1 is used. We adopt 
£’MeOH=o270 K as a representative value. This range 
is higher than published methanol physisorption ener- 
gies. For comparison, Sandford Allamandola give 
£m«oh = 4240K for pure methanol; this could not 
explain the W3(H 2 0) MeOH abundance if Td « 95 K. 
Similarly, the theoretical value of £ M eOH = 2060K re- 
quires T{\ as 38K - too low for W3(H 2 0) and predicts 
large methanol enhancements in cool cores which are 
not observed, e.g. at 55K in W3(IRS5) (Helmich and 
van Dishoeck 1997). 

Experimental measurements for E’hcooh are not 
available. Based on the trend of condensation tem- 
peratures of the pure substances (Nakagawa 1980) a 
relatively large value might be expected, much larger 
than the theoretical value of 2570K in Table 2. With 
a frequency factor of 2 x 10 l 2 s“ 1 , 100K dust, and 
an inferred gas:solid ratio of 100:1 in G34.3, equation 
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(20) yields £hcooh = 5867 K. 

Dimethyl ether will not participate in hydrogen 
bonding and the value of 2820 K given by Hasegawa 
k Herbst (1993) is probably appropriate. 

It should be emphasised that these estimates de- 
pend on the unknown ice abundances in these par- 
ticular cores. Note that the estimated £hcooh and 
EmcOH are for binding to the hypothetical OH-containing 
surface, silanol or organic residue, whereas once a wa- 
ter ice monolayer has formed Eh 7 o and the remaining 
E\ will be determined by their interaction with sur- 
face water molecules. 

5.2. Model 

The numerical model is that used by Charnley 
(1997a) except that the rate coefficients were updated 
from Millar et al. (1997) and the differential equa- 
tions for H2O , MeOH, HCOOH, Me20 and all other 
neutral molecules were modified in line with equa- 
tions (16) k (17), and including the binding energies 
of §5.1. Starting from an initial gas phase at 10K, 
comprising only of H2, He, H and their related ions, 
the dust temperature is raised to about 150K and the 
mantles are evaporated. It is assumed that no refrac- 
tory metals are present (see Charnley 1997b). The 
abundances (i.e. y(X)) of H2O , CO, MeOH, and 
HCOOH injected were taken from §5.1. The core has 
a density of na = 2 x 10 6 cm~ 3 , the visual extinction 
was taken to be 300 magnitudes, the cosmic ray ion- 
ization rate was taken to be 1.3 x 10"" 17 s _l . The 
effects of cosmic ray-induced photochemistry were ig- 
nored. 

Based on the binding energies of §5.1, it is as- 
sumed that the grains cool to where CH3OH and 
other molecules can recondense but that they remain 
hot enough that more volatile species, such as ethers, 
do not. To keep the model as simple as possible, a 
simple grain cooling law was adopted 

T d (<) = 150exp(— f/r COO |) (22) 

and calculate the core chemistry for various values of 
r coo i. In fact, r coo i is relatively well constrained by 
the requirements that it be less than the dynamical 
life-time of the core, about 2 x 10 5 years, but also be 
long enough for the CH3OCH3 abudnance to grow to 
its observed value. To produce CH3OCH3 /H2 abun- 
dances of ~ 10~ 8 for an injected CH3OH /H2 fraction 
of 10“ 6 places a lower bound on r coo | of around 5000- 
10,000 years (Charnley et al. 1995). In summary, the 


calculation proceeds with the dust instantaneously 
heated to 150K, as in previous hot core models, it is 
then allowed to cool to 100K on a time-scale of r coo \ 
where it is held fixed for 2 x 10 5 years. Calculations 
are presented for r COO i equal to 10 4 years, as well as 
for the case rcooi = 00 to account for the case of no 
recondensation, as assumed in all previous work. 

The models shown in Figure 1, where the calcu- 
lated abundances are presented as relative to molecu- 
lar hydrogen (i.e. 2y(X)), show that the recondensa- 
tion scenario can be realised in hot core models and 
that it has a dramatic effect on the core composition. 
It is suggested that this process is occurring to some 
extent in G34.3 and W3(H20). 

6. Conclusions & Discussion 

Early in the lives of embedded massive protostars, 
a slight reduction in the protostellar luminosity causes 
the local dust temperature to fall and permits recon- 
densation of molecules that were evaporated previ- 
ously when the dust was warmer. Such a scenario 
can explain the apparently anomalous abundances of 
methanol and its daughter molecules CH3OCH3 amd 
HCOOCH3 in G34.3 and W3(H 2 0), as well as the 
low abundances of HCOOH, and perhaps also EtOH, 
generally seen in hot cores. 

Hydrogen bonding must be important in determin- 
ing the binding energies of molecules containing the 
hydroxyl group and the appropriate values have been 
derived from the astronomical observations. This 
work suggests that it may be possible to determine 
molecule-grain binding energies for many species. If 
the physical conditions along particular sightline are 
well known, especially the dust temperature in the 
hot core component, then, for molecules produced on 
grains, the gasisolid abundance ratio will allow an es- 
timate of E\ to be made and depletions predicted for 
other cores. In practice one must be cautious since 
one is dealing with spatial gradients in density and 
temperature and so estimating values may not be as 
straightforward as presented here. 

One consequence of this theory is that the entire 
composition of interstellar ices will only be present in 
gas at 150K or more. The fact that many complex 
products of grain surface chemistry (e.g. alcohols) 
contain O-H and N-H bonds means that they will 
form hydrogen bonds with either a silanol-like surface, 
or organic residue, and consequently will have large 
binding energies. The fact that, relative to water, the 
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abundances of HCOOH, NHjCHO and HNCO in hot 
cores are significantly depleted with respect to those 
seen in the coma of Comet Hale-Bopp (see Crovisier 
1998) suggests that recondensation may be occurring 
in the former sources. 

Finally, molecular recondensation may also occur 
around low-mass protostars, either following mantle 
removal in shock waves or, in more evolved objects ,as 
a result of FU Ori outbursts. The latter case wil be 
discussed elsewhere and we note here that Avery and 
Chiou (1996) have already identified recondensation 
of CH 3 OH and SiO as a possible explanation for the 
peculiar chemical abundances observed in the LI 157 
bipolar outflow. 

Theoretical astrochemistry at NASA Ames is sup- 
ported by NASA’s Origins of Solar Systems and Exo- 
biology Programs. 
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Table 1 

OBSERVED ABUNDANCES IN HOT CORES 


Molecule 

Compact® 

Ridge 

W3(H 2 0) 6 

G34.3 c 

HjO d 

5000* 

770 

4300 

CH3OH 

94.0 

4.6 

10.0/3.4 

CH3OCH3 

4.7 

1.0 

1.0/- 

HCOOCH3 

7.8 

0.3 

1. 0/3.0 

HCOOH 

0.3 

- 

-/0.2 

T,«(K) 

70-150 e 

100 

75-190 

t»h (10®cm -3 ) 

> l e 

rv5 

1/ 


Entries are fractional abundances relative to H 2 divided 
by 10~®. Dust and gas are assumed to have equal temper- 
atures. 

* Sutton ct al. (1995); calculated using CO/H 2 =2.7 x 10~ 4 
(Helmich & van Dishoeck 1997) 

6 Helmich II van Dishoeck (1997) 

c For G34.3, two values are given: Mehringer & Snyder 


(1996) /Mac Donald et al. (1996) 
* Gensheimer et al. (1996) 


Blake et al. (1987) 

/ Millar et al. (1997b) for the 'compact core’. The ‘ultra- 
compact core 1 has ~ 300K 

1 Not detected specifically in the Compact Ridge; estimate 
based on HDO/H 2 = lCr 8 and HD0/H 2 0=2 x 10~ 4 


Table 2 

BINDING ENERGY DATA 


X 

X-SiOj° 

(K) 

X-H 2 O fc 

(K) 

x-x* 

(K) 

*x e 

(10“ 24 cm 3 ) 

CH30CH3 

2820 


- 

5.29 

h 2 o 

1860 

4815' 
5070 1 

4815* 

5070* 

1.45 

CH30H 

2060 

- 

4235 

3.29 

HCOOH 

2570 

- 

- 

4.42 

CO 

1210 

1740 

960 

1.98 

nh 3 

1110 

- 

3075 

2.26 

EtOH 

3470 

- 

- 

5.41 

C0 2 

2500 

2860 

2690 

2.63 


Binding energies denoted X-S are for adsorption of molecule 
X on surface S. a is the polarizability. 

tt Hasegawa Herbst (1993) and Allen Robinson (1977); 
theory. 

6 Sandford & Allamandola (1993); experiment. 

c Gonzalez et al. (1997) and Miller (1997); experiment and 

theory. 

1 unannealed ice. 

* annealed ice. 
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Fig. 1. — Chemical evolution in a hot core where 
molecules are recondensing on to cooling dust grains 
(full lines) and where no recondensation occurs (bro- 
ken lines). See text for details. 


This 2-colunm preprint was prepared with the A AS 
macros v4.0. 
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